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IS.  ABSTRACT 

A  method  is  developed  for  estimating  extreme  rainfall  rates  from  readily  avail¬ 
able  climatology  in  the  form  of  the  precipitation  index,  the  total  annual  rainfall  di¬ 
vided  by  the  number  of  days  with  rain.  The  method  is  expressed  as  an  equation  that 
estimates  the  rain  that  will  fall  during  any  selected  time  interval  (duration  of  rain¬ 
fall)  with  any  desired  frequency  of  occurrence  (return  period).  The  technique  as  pre¬ 
sented  is  applicable  to  the  geographic  regions  of  intense  precipitation  such  as  the 
moist  tropics  and  subtropics, 

To  provide  environmental  criteria  for  the  design  of  military  equipment  is  the 
objective  of  MIL-STD-210B.  The  equation  was  applied  to  this  problem  using  0.75  for 
the  precipitation  Index  as  typifying  the  areas  of  most  intense  rainfr.ll  for  calculation 
of  the  following  average  intensities  of  rainfall  in  inches  per  hour  that  would  be 
exceeded  with  a  probability  of  10  percent. 
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Abstract 

A  method  is  developed  for  estimating  extreme  rainfall  rates  from  readily 
available  climatology  in  the  form  of  the  precipitation  index,  the  total  annual  rain¬ 
fall  divided  by  the  number  of  days  with  rain.  The  method  is  expressed  as  an 
equation  that  estimates  the  rain  that  will  fall  during  any  selected  time  interval 
(duration  of  rainfall)  with  any  desired  frequency  of  occurrence  (return  period). 
The  technique  as  presented  is  applicable  to  the  geographic  regions  of  intense 
precipitation  such  as  the  moist  tropics  and  subtropics. 

To  provide  environmental  criteria  for  the  design  of  military  equipment  is 
the  objective  of  MI1.-STD-210B,  The  equation  was  applied  to  this  problem  using 
0.7  5  for  the  precipitation  index  as  typifying  the  areas  of  most  intense  rainfall  for 
calculation  of  the  following  average  intensities  of  rainfall  in  inches  per  hour  that 
would  be  exceeded  with  a  probability  of  10  percent. 
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Extremes  of  1,  12,  and  24  Hour  Rain 

for  MIL-STD-210B 


l.  INTRODI  CTIO.\ 

"Climatic  Extremes  for  Military  Equipment"  (MIL-STD-210B),  a  revision 
of  MIL-STD-210A,  will  present  logical  extremes  of  all  meteorological  elements 
that  could  have  an  impact  on  the  design  of  military  equipment,  IntenBe  rainfall 
can  affect  equipment  by  hampering  or  preventing  operation  while  the  rain  is  falling. 
Accumulation  of  unusual  amounts  can  cause  irreversible  damage,  so  that  the  equip¬ 
ment  is  worthless  for  further  operation  after  cessation  of  the  rain.  The  latter, 
termed  withstanding  extremes  for  design,  is  of  concern  in  this  study. 

To  withstand  rainfall,  equipment  must  be  able  to  survive  periods  of  intense 
rainfall  during  expected  dur  ations  of  exposure  (EDE)  in  the  field.  The  length  of 
the  periods  of  intense  rainfall  that  could  be  critical  and  the  EDE  will  vary  with 
the  equipment.  Hence  rainfall  intensity  data  must  be  provided  for  various  periods 
of  precipitation  and  EDEs,  Military  Standard  21013  specifies  duration  of  1,  12, 
and  24  hours  and  EDEs  shall  be  2,  5,  10,  or  25  years,  as  considered  appropriate 
for  each  item  of  equipment.  For  these  periods,  a  calculated  risk  of  failure  of 
10  percent  in  the  most  severe  geographical  area  for  each  climatic  element  is 
acceptable  for  "withstanding."  Data  on  extreme  annual  precipitation  are  available 
in  the  form  of  return  periods.  The  approximate  return  periods  that  will  give  a 
10  percent  chance  of  occurrence  within  the  specified  planned  life  spans  are: 

(Received  for  publication  17  May  1973) 


2 


5 


10 


25 


EDE  (years)  2 

Return  period  (years)  20  50  100  250 
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To  determine  values  of  rainfall  for  various  durations,  that  is,  1,  12,  and 
24  hours  for  specific  return  period,  extreme  probability  theory  (Gumbel,  1058; 
Gringorten,  19G3)  is  used.  Tabulation  of  annual  extremes  for  the.ee  durations  lor 
many  years  are  required.  Although  there  are  many  stations  throughout  the  world 
where  precipitation  is  measured,  there  are  relatively  few  where  amounts  are 
tabulated  and  published  for  time  periods  of  less  than  24  hours.  In  tropical 
regions,  where  the  most  intense  precipitation  is  to  be  expected  and  therefore  of 
greatest  military  interest,  published  data  on  annual  extremes  of  rainfall  are  even 
less  available  than  for  the  highly  industrialized  mid -latitudes.  In  order  to  provide 
climatological  information  on  a  world-wide  basis,  it  is  necessary  to  relate  data 
available  on  extreme  rainfall  for  mid-latitudes  to  the  usual  meteorological  obser¬ 
vations  for  which  a  world-wide  climatology  is  available. 

Rainfall  data  suitable  for  establishing  such  a  relationship  are  published  in 
l 'SUB  Technical  Paper  No.  25  (1955)  which  presents  data  for  200  stations  in  the 
contiguous  United  States  plus  1  station  each  in  Alaska,  Hawaii,  and  Puerto  Rico, 
Data  are  presented  as  rainfall  intensity-duration  curves  for  return  periods  of 
2,  5,  10,  25,  50,  and  L00  years  by  station,  Duration  ranges  from  5  minutes  to 
24  hours.  The  curves  are  based  on  the  maximum  rainfall  for  each  year  for  each 
duration.  Curves  for  the  return  periods  are  spaced  according  to  the  Fisher -Tippet 
type  I  (Gumbel)  distribution.  As  an  example  of  the  form  of  data  presentation,  the 
set  of  curves  for  Pensacola,  Florida  is  reproduced  in  Figure  1,  This  is  the 
station  with  the  most  intense  rainfall  of  the  203  available.  Rainfall  amounts  were 
extracted  for  all  of  the  stations  for  a  range  of  dur  ations  nud  return  periods.  Tire 
durations  and  periods  used  in  the  analysis  are  shown  in  Table  1,  which  also  con¬ 
tains  the  average  r  ainfall  rates  for  Pensacola  as  extracted  from  Figure  1, 

Similar  intensity -duration  curves  are  available  for  only  a  limited  number  of 
locations  outside  of  the  United  States,  a  few  in  areas  of  high  rainfall.  Included 
in  this  study  are:  Hong  Kong  (Cheng  and  Kwok,  1966);  Bombay  (Patel  and  Vaniari, 
1969);  Nagasaki  and  Yokohama,  Japan;  Naha,  Okinawa;  Manila,  San  Fernando  and 
Olangapo,  P.I,  (Paulhus,  1964). 
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The  objective  of  the  analysis  is  to  relate  extremes  of  I,  1 ami  24  hour  rain¬ 
falls  to  elimatolugical  data  that  are  widely  availtible.  The  simplest  hypothesis  is 
that  the  greatest  extremes  will  occur  in  areas  of  greatest  total  precipitation. 

Indeed  there  is  a  correlation  between  annual  pi  oeipitation  total  am  the  intensity  o' 
extreme  rainfall  that  is  statistically  significant.  The  relationship  is  of  limited 
practical  value,  however,  as  the  cor  relation  is  not  high,  and  the  standard  error  of 
estimate  of  the  regression  equation  is  fairly  large,  A  more  sophisticated  hypothesis 
is  one  t hat  would  take  account  of  Ihe  distritiul ion  of  the  average  dailv  iatrnsitv  of 
rainfall.  This  may  be  achieved  by  using  the  uerage  pro  ipd.Unm  on  a  rain1,  dav, 
obtained  by  dividing  the  annual  precipitation  by  the  number  of  days  on  which 
precipitation  occurred.  This  parameter  has  been  designated  as  t tic  precipitation 
index  (I)  and  the  primary  independent  variable  in  the  relationships  that  have  been 
developed.  It  is  fairly  high  in  correlation  with  extreme  intensities,  and  t tic 
relationship  has  a  standard  error  that,  although  larger  than  desirable,  is  small 
enough  to  have  practical  value  for  estimated  intensities. 

Other  predictors  were  tried,  alone  ami  in  multiple  regressions.  Only  those 
involving  the  frequency  of  thunderstorms  and  temperature  yielded  any  significant 
results.  Several  variables  were  created  io  serve  as  indices  of  thunderstorm 
activity,  for  example,  number  of  days  with,  thunderstorms/number  of  days  with 
rain.  Some  variables  improved  the  intensity  estimates  considerably  and  were 
statistically  significant.  However,  these  were  abandoned  with  great  reluctance 
owing  to  an  anticipated  lack  of  availability  of  thunderstorm  climatology  and  uniformity 
of  observations,  It  was  believed  that  there  would  be  areas  in  wh:ch  stations  pro¬ 
viding  thunderstorm  climatology  would  be  scarce  arid  that  wlici  c  data  we  re  avail¬ 
able,  there  would  be  differences  between  regions  in  the  climatology  presented. 
Specifically,  the  ecunt  of  number  of  days  w  ith  thunderstorms  was  suspect,  i  he 
authors  believed  it  quite  likely  tha,  this  statistic  would  not  tiavo  tiie  same  meaning 
in  all  regions  of  the  world. 

The  teiiipvruture  parameter  found  by  trial  to  be  most  useful  was  an  average 
measure  of  the  annual  temperature  range:  the  difference  in  mean  monthly  tempera¬ 
ture  in  b'ah  "enheit  degrees  bet  ween  the  wannest  and  coldest  months.  Tins  param¬ 
eter  UIT)  was  taker  as  the  secondary  independent  variable, 

Regression  equations  wore  established  for  each  return  period  and  duration 
based  on  data  front  all  stations.  A  few  of  these  equations  were  ehiekrt! 
against  several  of  the  stations  with  most  intense  precipitation.  There  seemed  to 
be  some  slight  tendency  to  underestimate  these  extreme  values.  Hence  it  was 
decided  to  select  a  smaller  sample  of  stations  with  high  ra.es  oT  precipitation.  In 
the  United  States,  these  stations  were  located  near  the  coast,  from  South  Caroline 


*** 


southward  on  the  Atlantic  coast  and  all  along  the  Gulf  coast  from  Brownsville, 

Texas  to  Key  West,  Florida.  All  previously  listed  stations  outside  of  the  United 
States  were  included  in  this  smaller  sample,  consisting  of  27  stations. 

Regression  equations  were  established  based  on  this  sample.  The  correlations 
were  lower  than  for  tne  regressions  based  on  all  stations  but  so  was  the  standard 
error  of  estimate.  For  this  email  sample,  multiple  regression  was  of  no  value; 
the  contribution  to  explained  vat  iance  of  the  variable  dT  was  not  statistically  sig¬ 
nificant,  nor  was  the  standard  error  of  estimate  improved  appreciably.  Regressions 
utilizing  all  stations  were  improved  by  the  addition  of  dT  as  an  independent  variable 
in  all  but  'he  12  and  24  hour  durations.  The  standaro  error  of  rates  estimated  from 
these  regressions  was,  however,  larger  than  for  estimates  made  from  regressions 
on  (I)  only,  derived  from  the  sample  of  27  stations.  Since  these  27  stations  have 
the  most  intense  precipitation  and  provide  the  greater  precision  of  estimate,  the 
model  to  be  developed  will  be  based  on  them. 

The  estimating  equation  to  be  used  is 

R  =  A  +  BI  ,  <n 

where  R  is  the  rair  rate  in  inches  per  hour  and  I  is  the  precipitation  index  in 
inches  per  day  of  rainfall  >  0.01  inch.  Since  there  are  8  durations  (D)  and  4  return 
periods  (P),  the  equation  can  be  rewritten  with  subscripts  to  denote  the  particular 
duration  and  return  period: 

RUP  =  ADP  +  BDP!  ’  (la) 

The  32  regression  equations  were  determined  by  least  squares  and  the  co¬ 
efficients  are  shown  in  Figures  2  and  3  on  a  semilogarithmic  scale.  They  appear 
to  increase  linearly  with  the  logarithm  of  P,  the  return  period.  This  relationship 
is  expressed  by 


where  is  the  coefficient,  either  A  or  B;  for  example,  A^  -  +  b^glnP  . 

This  model  was  fitted  oy  least  squares  to  the  16  sets  of  values  (2  coefficients 
X  8  durations)  with  all  fits  being  highly  significant  statistically.  The  lowest  cor¬ 
relation  obtained  was  0,  967  whic'.  provided  an  F  ratio  that  was  significant  at  the 
3  percent  level.  The  values  of  the  coefficients  and  bj^  are  shown  in  Figure  4 
on  a  log-log  scale.  There  appears  to  bo  a  quadratic  relationship  between  the  value 
of  a  coefficient  and  D,  the  duration.  This  relationship  is  expressed  by 


Figure  2.  Intercept  of  Regression 
Lines  Fitted  to  27  Rain  Rates  for  Dura¬ 
tions  Specified 


Figure  3.  Slope  of  Regression  Lines 
Fitted  to  27  Rain  Rates  for  Durations 
Specified 
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Figure  4.  Coefficients  of  C.  -  a.  -*  b.  In  F 
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In  kb  -  o  +  ;ij  InD  +  y.  ln2D  ,  (3) 

where  kb  is  the  coefficient  or  in  Eq,  u,,  ■  example,  In  ^  + 

3a  In  D  +  yA  ln2D  . 

This  model  was  fitted  to  the  4  sets  of  coefficients  yielding  the  following  four 


equations,  all  with  c<>rrc;r,.ions  greater  than  0.99. 

In  A  -  1 . 33  1 V  *  0.  22135  In  D  -  0.  13880  ln2D  (4a) 

a 

In  A,  -0.4 11243  J  0.  33 f.52  In  D  -  0.094112  ln2D  (4b) 

u 

In  B  --  0.  58770  +  0.  15912  In  D  -  0.05524  ln2D  (4c) 

a 

In  Ub  0.  81299  -  0.  (12919  In  D  -  0.  02514  ln2D  (4d) 

Equations  (1)  and  (2)  can  be  combined  into  a  general  equation, 

R  =  (A  +  A.  InP)  4-  (II  +  B,  InP)  I  ,  (5) 

cab  a  b 


and  the  values  of  the  coefficients  A  ,  A,  ,  B  and  Bu  determined  from  Eqs,  (4), 

a  b  a  b 

Here  R  ,  is  used  to  designate  t he  computed  or  estimated  value  of  precipitation  ns 
distinguished  from  It,  the  observed  value  that  was  the  input  to  the  analysis. 
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1.1  l.lfcrl  of  Nn»olhiii|!  Hrprr'Moo  t.orfficii'HlN 

Ordinarily  with  least-squares  regression,  it  is  possible  to  make  a  statement 
of  ihe  precision  of  the  estimate  yielded  by  the  equation  by  presenting  the  standard 
error  of  estimate.  This  can  be  done  for  the  initial  stage  of  32  equations  (1), 

When  equations  (2)  are  determined,  the  standard  error  obtained  applies  to  the 
estimated  values  of  coefficients  in  equations  (1),  not  to  the  values  of  R  estimated 
from  those  equations. 

In  turn,  the  use  of  Eq.  (3)  to  obtain  estimates  of  coefficients  in  Eq.  (2)  does 
not  give  a  measure  of  the  precision  of  estimate  for  rainfall. 

The  only  way  to  obtain  a  statement  of  the  precision  of  estimating  R  from  the 
smoothed  regression  coefficients  is  by  calculating  the  actual  individual  differences, 
R-Rc,  between  observed  and  estimated  values  for  each  station,  duration,  and  return 
period.  These  differences  can  be  accumulated  and  compared  to  the  standard  error 


obtained  from  the  basic  regressions  (Eq.  1).  For  purposes  of  this  comparison, 
the  statistic  to  accumulate  would  be  the  root-mean-square 


f  2  "]  2 

RMS  =  |<R-Rcr/NJ 

In  error  theory  this  is  a  statement  of  precision,  a  statement  of  variation  about  a 
zero  mean  with  positive  and  negative  deviations  being  eaual  in  number  and  mag¬ 
nitude.  The  standard  error  of  estimate  in  regression  is  an  analogous  statement 
of  variation  about  the  regression  curve,  with  positive  and  negative  deviations 
being  equal  in  number  and  magnitude.  In  this  case,  however,  the  deviations  are 
about  a  curve  that  is  no  longer  the  least-squai  ta  regression  so  that  positive  and 
negative  deviations  are  no  longer  equal.  In  fact,  there  can  be  a  net  bias  and  the 
mean  of  the  deviations  can  be  other  than  zero.  One  could  compensate  for  this  by 
calculating  the  standard  deviation  o,  a  measure  of  the  dispersion  of  the  deviations 
about  their  mean.  However,  this  is  not  a  measure  of  the  precision  of  the  estimate 
In  order  to  accomplish  an  evaluation,  both  the  RMS  and  a  values  were  cal¬ 
culated,  along  with  the  mean  deviation  for  each  of  the  3  2  combinations  of  duration 
and  return  period.  Those  are  presented  along  with  the  standard  error  of  estimate 
from  the  b.»sic  regressions  {Eq.  i)  in  Table  2.  Fortunately,  the  mean  deviations 


Table  2,  Measures  of  Precision  of  Estimate  * 


\  p 

Standard 

Error  (in. 

/hr) 

RMS  (in.  /hr) 

D  \ 

2 

10 

50 

100 

2 

10 

50 

100 

\ 

(years) 

(years) 

5  min 

0.  533  7 

0.7515 

1.  1711 

1.3170 

0. 5390 

0.  7532 

1. 1720 

1 . 3  1 7  G 

1  5  min 

0.  3939 

0.  5467 

0.  7655 

0.  7741 

0,  4052 

0 . b‘ 7  0 

0.  7G7i 

0.  7771 

30  min 

0.  273  7 

0.4415 

0.  5785 

0.  5794 

0.  2771 

0.4443 

0. 5703 

0.  5839 

1  hr 

0.  1653 

0.3032 

0.4144 

0.4385 

0.  1661 

0.3113 

0.4150 

0.4404 

3  hr 

0.  13  59 

0.  1857 

0.  2680 

0.3188 

0.  1378 

0.  194  5 

0.  2  706 

0.  3219 

(i  hr 

0.  1018 

0.  1458 

0.  2005 

0. 2268 

0.  1026 

0.  1528 

0.  2042 

0.  2304 

12  hr 

0.  073  9 

0.  1072 

0.  1567 

0.  1810 

0.  0740 

0.  1079 

0.  1568 

0.  1810 

24  hr 

0.  0493 

0.0734 

0.  1070 

0.  1194 

0.  0505 

0.0736 

0.  1078 

0.  1202 

ct  (in 

/hr) 

Mean  (in.  /  hr) 

5  min 

0.  5340 

0,7517 

1.  171  1 

1.3173 

-0.  0732 

0.0475 

-0.0459 

0.  0278 

1  5  min 

0.4004 

0.  5469 

0.  7664 

0.  7741 

-0,  0625 

•0.0088 

0.  0332 

-0. 0679 

30  min 

0.  273  9 

0.4432 

0.  5786 

0.  5796 

-0.0420 

0,0312 

-0. 0279 

-0,  0704 

1  hr 

0.  1656 

0.3055 

0.4149 

0.4394 

0,0133 

0.  0601 

-0.0097 

-0. 0291 

3  hr 

0.  13  67 

0.  1858 

0.  2680 

0.3204 

0.  0175 

0.0574 

0.0374 

0.  0305 

6  hr 

0.  1024 

0.  1466 

0. 2008 

0.  2270 

0.  007  5 

0.0431 

0.037a 

0.  0389 

12  hr 

0.0739 

0. 1073 

0.  1568 

0.  1810 

-0.0047 

0.012n 

-0,0012 

0.  0012 

24  hr 

0.0493 

0.0735 

0.  1071 

0.  1195 

-0.0108 

-0.  0016 

-0,0129 

-0.  0135 

*  Standard  Error  of  basic  regressions  [  Eq.  (1)1;  other  measures  calculated 
from  deviations  (R-Pc)  between  actual  values  and  estimates  from  equations  with 
smoothed  coefficients  for  27  station  subset. 
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from  the  equations  with  smoothed  coefficients  are  smell  so  that  there  Is  little  dif¬ 
ference  between  the  values  of  a  .  nd  RMS.  In  fact,  the  RMS  is  only  slightly  larger 
than  the  standard  error  of  estimate  of  the  basic  regressions.  Hence  it  is  concluded 
that  the  calculation  of  orecipltation  rates  based  on  Eqs.  (4)  ann  (5)  is  a  valid  pro¬ 
cedure. 

1.2  Kffect  of  I. imi line  Sample  Size 

If  the  complete  set  of  211  stations  is  used  to  determine  the  regression  equations, 
the  variable  dT  contributes  significantly  to  the  explained  variance.  T  he  subset  of 
27  stations  indicated  that  dT  did  not  contribute  sig”i(. candy.  Thus  there  are  two 
plausible  approaches  to  using  the  information  from  the  full  data  set  in  estimating 
extreme  precipitation  rates:  A  simple  relationship  to  precipitation  index  (1)  or  a 
multiple  relationship  to  both  I  and  dT.  Both  sets  of  equations  were  used  and  the 
deviations  (n-Rc)  between  the  observed  and  computed  values  for  the  27  stations 
subset  were  compiled  into  RMS  statistics  which  are  presented  in  Table  3.  Nearly 
all  of  the  values  in  Table  3  are  appreciably  larger  than  the  RMS  values  based  on 
equations  derived  from  the  subset  as  shown  in  Table  2.  Hence  it  is  concluded  that 
limiting  the  analysis  to  the  27  stations  with  extreme  precipitation  rates  is  profit¬ 
able. 


Table  3.  RMS  Deviation  (R-Rc)  Between  Actual  and  Estimated  Values  for  27  Station 
Subset  Based  on  Equation0'  With  Smoothed  Coefficients  Derived  From  All  Sations 


I _  _  . 

Nf 

1  only 

I  and  dT 

'D 

1 _ X 

2 

10 

50 

100 

2 

10 

50 

100 

1  5  min 

1.3227 

1.  1266 

1.3608 

1. 5016 

1.2770 

1.  5691 

1.9147 

2.0834 

1  5  min 

0. 8563 

0.9679 

1.  1051 

1.  1458 

1.0474 

1.3991 

1.7540 

1.8541 

|  30  min 

0.  6340 

0.8330 

0. 9471 

0.  5779 

0.7860 

1.  1575 

1.4383 

1,  5231 

|  1  hr 

0.4232 

0.  6541 

0.7909 

0.  8552 

0.  5432 

0.  8593 

1.0753 

1.  1706 

3  hr 

0.  2024 

0.3822 

0. 5293 

0, 5862 

0.  2410 

0.4426 

0. 6082 

0.  6712 

6  hr 

0.  1292 

0.2484 

0.3583 

0.4084 

0. 1354 

0.  2703 

0.3951 

0.4520 

12  hr 

0.0893 

0.  1632 

0.2346 

0.2702 

0.0824 

0.  1600 

0.  2401 

0.  2794 

24  hr 

0.0643 

0.  1166 

0.  1608 

0.  1813 

0.0599 

0.0921 

0. 1227 

0.  1363 

\ppijcvnos 


5. 1  Precision  of  Estimate 

Whatever  the  application,  a  knowledge  of  the  precision  of  the  estimate  may  be 
desired.  This  may  be  obtained  from  the  RMS  section  of  Table  2.  The  usual  inter¬ 
pretation  of  the  RMS  statistic  is  based  on  the  assumption  of  a  normal  frequency  dis¬ 
tribution  with  zero  mean.  The  frequency  distributions  of  the  deviations  (R-Rc) 
were  checked  for  skewness  and  kurtosis.  A  few  of  the  32  sets  were  nearly  normal 
but  most  departed  from  normal  but  not  by  great  amounts.  There  was  no  pattern 
for  the  departures.  It  is  probable  that  no  great  error  would  be  made  in  accepting 
the  usual  interpretation  of  the  statistic;  that  is,  that  the  true  value  of  R  lies  in  the 
range  R£  ±  RMS  with  a  probability  of  63  percent. 

When  the  regression  equations  are  extended  beyond  the  range  of  the  data, 
either  in  duration  or  return  period,  Table  2  will  not  provide  an  estimate  o !  the  RMS 
deviation.  The  values  listed  in  Table  2  show  an  orderly  progression  as  did  the 
coefficients  in  the  regression  equations.  The  RMS  values  are  linearly  related  to 
the  logarithms  of  the  return  period  as  can  be  seen  in  Figure  5.  The  relationship 
is  expressed  by 

RMS  =  A  +•  B  In  P  (6) 


Figure  5.  RMS  Deviation  <R-RC) 
Between  Actual  end  Computed  Precipita 
tioii  Rates  for  Sample  of  27  Stations 


to 

« 


1 

o 

o 


Figure  6.  Coefficients  of  RMS  =  A  +  B  In  P 


which  was  fitted  by  least  squares  and  the  behavior  of  the  coefficients  A  and  B 
examined.  A  linear  relation  appears  in  Figure  6  where  A  and  B  are  plotted 
against  duration  on  a  log-log  scale.  Least- squares  regression  yielded  the  follow¬ 
ing  expressions  for  A  and  B  in  Eq.  6: 

In  A  =  -0. 14536  -  0.42516 InD  (7a) 

In  B  =  -1.06942  -  0.39679 InD  (7b) 


•1.2  Military  Standard  2 1  OH 

To  obtain  values  of  R  for  M1L-STD-2 1CB  using  Eqs,  4  and  5,  a  value  must 
be  assigned  to  the  precipitation  index  (I).  This  should  be  representative  of  the 
most  severe  geographical  area.  The  areas  with  most  intense  rainfall  are  found 
in  Africa  on  the  Gold  Coast,  in  South  America  in  Columbia  and  in  the  Philippines, 
Indonesia,  South  East  Asia  and  India,  The  annua),  precipitation  index  was  cal¬ 
culated  for  a  number  of  stations  in  the  Asian  region  which  has  the  largest  area 
of  intense  rainfall,  A  few  stations  have  an  index  greater  than  1  inch  per  rain 
day  of  0,01  inch  or  more,  but  these  are  relatively  rare  and  isolated.  Most  of 
Southeast  A3ia,  India,  Indonesia  and  the  Philippines  have  an  index  greater  than 
0,  5  in.  per  day.  A  number  of  stations  in  Burma  and  Southeast  Asia  yielded 
index  values  above  0.75  in,  per  day  and  formed  coherent  regional  patterns.  This 
index  value  appeared  to  be  representative  of  area3  of  greatest  rainfall  intensity 
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without  being  of  such  unuaual  occurrence  as  to  produce  an  unrealistically  stringent 
design  criterion  that  would  cause  excessive  overdesign  of  equipment.  Hence 
I  =  0.75  was  selected  for  input  into  Eqs.  4  and  5  for  calculating  values  of  Intense 
rainfall  for  MIL-STD-210B. 

For  the  purposes  of  MIL-STD-210B  the  durations  of  interest  are  1.  12  and 
24  hours  and  the  return  periods,  as  stated  previously  are  20,  50,  100  and  250  year 
The  coefficients  obtained  from  Eqs.  4  and  2  are  given  in  Table  4. 


Table  4.  Coefficients  for  Calculating  Rain  Rates  for  M1L-STD-210B 


Coefficient 

Duration  (hours) 

1 

12 

24 

A(a) 

0.  9132 

0.03  98 

0.0122 

A(b) 

0. 5113 

0.0959 

0.0488 

B(a) 

1.3678 

0.4693 

0.3083 

Bib) 

0.2514 

0.0966 

0.0778 

20  yr  A 

2.4450 

0.3272 

0.  1585 

B 

2.  1209 

0.7588 

0.  5413 

50  yr  A 

2.9135 

0.4151 

0.2033 

B 

2.3512 

0.8473 

0.  6125 

100  yr  A 

3.2679 

0.4815 

0.  2371 

B 

2.  5255 

0.9143 

0.6664 

250  yr  A 

3.7364 

0.  5694 

0.2819 

B 

2.7559 

1.0028 

0.7377 

5.3  Definition  of  I’rci'ipitulion  Index 

The  precipitation  index  (I)  has  been  defined  as  the  total  annual  precipitation 
in  inches  divided  by  the  number  of  days  with  0.  01  inch  of  precipitation.  This 
was  chosen  because  the  climatological  records  for  almost  all  of  the  stations 
used  in  this  study  contained  the  parameters  in  these  units.  There  is  no  problem 
in  converting  total  annual  rainfall  from  millimeters  to  inches.  Unfortunately, 
not  all  nations  define  a  day  with  precipitation  in  the  same  manner.  Some 
count  days  with  l  mm  or  0,  1  mm  or  1  /  10th  inch.  When  these  are  encountered 
it  is  necessary  to  convert  to  the  number  of  days  with  0,01  inch  in  order  to  cal¬ 
culate  the  index  for  use  in  Eq.  (5).  Conrad  and  Poliak  (1950)  present  a  technique 
for  doing  this  by  calculating  the  probabilities  p0,  p^,  p^,  ...,  p.  of  days  with 
0,  1,  2,  . . .  i  units  of  rainfall.  This  is  described  by 


Pi 


_  h  +  (i-1)  d 

pi-i  rrr+ar 


and 
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P°  (l+d)h/d  ' 

where  h  =  RR/N,  RR  Is  the  total  rain/all  for  the  period  In  the  units  specified  and 
N  is  the  total  number  of  days  involved.  This  would  be  365.  25  when  dealing  with  a 
yearly  average  record.  Thus  h  is  a  fictltioua  average  rainfall  for  one  day.  The 
variable  d  represents  the  degree  of  dependence  of  one  rainfall  event  on  another. 

It  must  be  determined  from  the  ri  -d  which  will  give  the  count  nj  of  the  number  of 
days  with  rainfall  equal  to  or  gr  .han  i  units.  Then  (N-n^/N  is  the  probability 

of  less  than  i  units  of  rain  occum .  This  probability  is  also  given  by 

2  Pj.  Trial  values  of  d  are  assumed,  £p^  is  calculated  and  compared  to  (N-n^l/N 

until  a  value  of  d  is  found  that  equates  the  two  estimates.  Note  that  the  units  in¬ 
volved  are  not  millimeters  or  Inches.  When  converting  from  days  with  0.  1  in.  or 
more  of  rainfall,  one  finds  that  the  unit  is  0.01  in.,  i  =  10  and  an  annual  total  of 
23.41  inches  is  2341. 


6.  SUMMARY 
6.1  Method 

A  method  of  estimating  extreme  rainfall  amounts  likely  to  occur  during  a 
specified  time  interval  D  with  a  return  period  P  from  widely  available  climatological 
data  has  been  developed.  The  climatological  input  is  the  precipitation  index  (I) 
obtained  by  dividing  the  total  annual  rainfall  in  inches  by  the  number  of  days  with 
0.01  Inch  or  more  of  rain.  The  method  is  described  by  the  .ollowing  equations 
where  R  is  the  estimated  rainfall: 

Rc  =  A  +  B  I 

A  -  A  +  A.  In  P 
a  b 

B  =  B  +  B,  In  P 
a  b 

In  A  =  1.33123  +  0.  22135  In  D  -  0,  13889  ln2D 
a 

In  Afa  =  -0. 46243  +  0. 33652  In  D  -  0. 09462  ln2D 

In  B&  =  0,  58770  +  0.  15912  In  D  -  0.  05524  ln2D 

In  Bfe  =  0.81299  -  0.  62919  In  D  +  0.  025)4  ln^D 

The  standard  error  of  estimate  (RMS)  of  R  from  the  foregoing  equations 
may  be  approximated  by 


u 


RMS  =  a  ♦  0  In  P 
where 

lna  ■  -0.  14536  -  0.42516  In  D 
In  0  •  -1.  06842  -  0.  38679  In  D 

The  model  Is  restricted  in  applicability  to  areas  of  intense  rainfall  such  as 
the  tropics. 

6.2  Application 

For  M1L-STD-210B  a  precipitation  index  value  of  0.75  was  used  to  calculate 
the  following  average  intensities  of  rainfall  in  inches  per  hour  that  would  be 
exceeded  with  a  probability  of  10  percent: 


Duration 

(hours) 


Estimated  Duration  of  Exposure  (years) 
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Station 

Index 

Station 

Apalachicola,  Fla. 

0.530 

Nagasaki,  Japan 

Austin,  Tex. 

0.402 

Naha,  Okinawa 

Bombay,  India 

0.711 

New  Orleans,  La. 

Brownsville,  Tex. 

0.377 

Olangapo,  P.  I. 

Charleston,  S.  C. 

0.427 

Pensacola,  Fla. 

Del  Rio,  Tex. 

0.300 

Port  Arthur,  Tex. 

Galveston,  Tex. 

0.498 

San  Antonio,  Tex. 

Hong  Kong 

0.597 

San  Fernando,  P.] 

Houston,  Tex. 

0.450 

San  Juan,  P.R. 

Jacksonville,  Fla. 

0.464 

Savannah,  Ga, 

Key  West,  Fla. 

0.357 

Tamna,  r'la. 

Manila,  P.I. 

0.  560 

Thomasville,  Ga. 

Miami,  Fla. 

0,474 

Yokohama,  Japan 

Mobile  Ala. 

0.  549 

Appendix 

Lilt  of  Station*  In  Samplo 

Index 
0.463 
0.499 
0.  534 
1,222 
0.  566 
0.  510 
0.362 
0.853 
0.307 
0.453 
0.460 
0.440 
0.408 


